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We investigate the possibility of having horizons inside various classical field 
configurations. Using the implicit function theorem, we show that models satisfying 
a certain set of criteria allow for (at least) small horizons within extended matter 
fields. Gauge and global monopoles and Skyrmions satisfy these criteria. Q-balls 
and Boson stars are examples which do not and can be shown not to allow for 
horizons. In examples that do allow for horizons, we show how standard 'no hair' 
arguments are avoided. 
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1. Introduction 



Black holes are intriguing objects and worth studying in all their possible 
varieties. In this paper we will study the possibility of having black holes inside 
various classical field configurations. Examples we consider include gauge and 
global monopoles, Skyrmions, Q-balls [1], and Boson stars [2]. 

Besides the basic search for black hole solutions, there are a number of physi- 
cally motivated questions one can ask in this context. For instance, what happens 
when you drop such an object into a Schwarschild black hole? For a gauge or global 
monopole the result should be a black hole with the appropriate kind of hair, since 
these both involve non-trivial behavior of the fields at infinity. But in the case of a 
Skyrmion, one might think that the only possibility would be its vanishing without 
a trace. Our results show that there is another possibility, at least for horizons 
very small compared to the Skyrmion radius^. In the case of gauge monopoles, 
Lee et.al. [7] have argued that, besides the Reissner-Nordstrom type solutions [8], 
there also exist, for sufficiently small horizon radius, solutions in which the Higgs 
field and gauge field behave more like an extended monopole outside the horizon. 
Our results confirm their arguments, and show that global monopoles can also have 
horizons inside them. 

Horizons inside extended field configurations may also be relevant in the late 
stages of black hole evaporation by Hawking radiation. Lee et.al. [9] have shown 
that extreme, magnetic Reissner-Nordstrom type black holes are unstable in a the- 
ory with extended monopole solutions. They conjecture that the extended solution 
discussed above is stable and that evaporation of the black hole proceeds through 
this configuration, leaving a non-singular magnetic monopole as the end state. 
Perhaps a Schwarschild black hole in a Skyrmion theory, for example, similarly 
becomes unstable (or metastable) when its radius is less than the characteristic 



f numerical results on extended Skyrmion fields around a black hole are given in references 
[3,4,5,6]. 
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Skyrmion radius. The evaporation process may then leave behind other stable 
remnants. 

Finally, in the literature Q-stars (large Q-balls) [10,11] and Boson stars (see 
[12] and references therein), as well as strange matter [13] and other types of non- 
topological solitons, are discussed as candidates for compact astrophysical objects. 
We can ask what the possible final collapsed states of such matter are. 

2. Existence of Solutions with Horizons 

We will be looking for static, spherically symmetric solutions to Einstein's 
equation, which have nonsingular, nontrivial matter fields outside a horizon. The 
form of the metric will be taken to be 

ds 2 = -B(r)dt 2 + A{r)dr 2 + r 2 dtt 2 (2.1) 

It is often convenient to define the function m(r) by 

1 2Gm(r) 
A(r) r 

A horizon occurs at coordinate rjj if 

2Gm{r H ) = r H . (2.3) 

When a horizon is present, one also expects that m = m(0) ^ 0, so that the 
metric is not well behaved at the origin. This is like having a seed mass at the 
origin. 

Let us agree to call a star a configuration of matter fields <fi (not necessarily 
scalar) such that the stress-energy is static, spherically symmetric, and localized. 
Suppose a particlar matter field theory has star type solutions, without gravity. 
There is some force balance, without gravity, which keeps the field configuration 
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(2.2) 



from either collapsing to a point or expanding to infinity. One might expect that 
weakly gravitating solutions would then exist, and that even placing a small seed 
mass inside the star, wouldn't disturb the balance too much. This can be made 
more precise by considering the Oppenheimer-Volkoff (OV) equation of hydrostatic 
equilibrium, which states (in the case when the three principal pressures are not 
necessarily the same) 

dpf „(m(r) + Anr 3 pf) 2. , 

dr ^ r(r-2G,n) (P + Pf) + r {P t- Pf) ^ 

In the absence of gravity, only the second term on the right hand side is present 
and for weak gravity, this term may still dominate. However, from the first term 
in (2.4), we see that at a horizon the sum of the radial pressure and the energy 
density must vanish. In a normal, burning star, both these quantities are positive 
and a horizon is not possible. On the other hand for many field theories, it happens 
quite naturally that (p + Pf)\ rH = 0. 

Our main result will be to show that given a matter theory which (1) has star 
type solutions without gravity and (2) satisfies (p + Pf)\ rH = "automatically" 
(in a sense defined below), then there exist star solutions when the matter theory 
is coupled weakly to gravity, and there also exist solutions with horizons inside. 

More precisely, the non-gravitating matter theory is described by a Lagrangian 
L m . A star solution is found by evaluating the action on field configurations 
consistent with a particular static, spherically symmetric ansatz. The Lagrangian 
restricted to this class of fields will be written L m (0). We will assume that — L m ((f>) 
is positive definite. When the matter theory is coupled to gravity, we will assume 
that the sum of the energy density and radial pressure is given by 

(pr + p) = jK(<i>), (2.5) 

Where K is a functional of the matter fields only. Then there exist regular star 
type solutions to the Einstein equation, and there also exist star type solutions 
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horizons, which have nontrivial, nonsingular matter fields outside the horizon, for 
G and r# sufficiently small. The argument, as follows, is an application of the 
implicit function theorem. 

First define new gravitational variables, 

e x = ^ and e y = VAB. (2.6) 
The action for fields outside a horizon is then taken to be S = Se + S rn , with 

oo 

Se(x,v) = J dry'((r -r H )e y -re x ) 

(2-7) 

S m ((p n ,x,y) = J drr 2 e y L m 

Se differs from the usual Einstein action by a boundary term, which has been 
chosen so that varying Se imposes the correct boundary condition at the horizon 
(see reference [14]). Varying the action with respect to x and y gives the equations 
of motion 

y > = ^ Gre y^ 5 llIL (2.8) 
ox 

Ur{e y - e*)) = -rV(e^^ + L m + ^) (2.9) 
ar ox oy 

and the boundary condition 



Note that from the definition of the stress tensor —2^- = p? + p. Equations (2.8) 
and (2.9) can be used to solve for the gravitational fields x and y in terms of the 
matter fields alone if and only if + p — j^K ((/)), where K is a function of the 
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matter fields alone. This was one of our assumptions. This is equivalent to the 
matter lagrangian having the form 

L rn (<f>) = -jK(<f>)-U(<f>,AB) (2.11) 

We can then define a positive definite functional of the fields, E(<j>, G, r#), by 

oo 

E{<j>\ G, r H ) = -S = J dre y (r(r - r H )K + r 2 U) (2.12) 
In (2.12), y(r) is given in terms of the matter fields by 

00 

y (r) = -8nG J dr'r'K(<p) (2.13) 

r 

Note that for a given configuration of the fields 0, E(<j>, G, r#) is a continuous, 
differentiable function of G and rjj- 

We assume that for G = rjj = 0, the functional E((/>, 0, 0) has a minimum O . 
This is our non-gravitating star. For G and r# nonzero, we seek solutions to 

F((&G,r H )) = ^ = 0, (2.14), 

which by construction will satisfy the equation of motion with the correct boundary 
conditions. By assumption F(<f) ; 0, 0) = 0. The implicit function theorem for 
Banach spaces* [15] can then be used to show that for G and rji sufficiently close 
to zero, there exist functions 0(G,r/J satisfying (2.14), such that 0(0,0) = O . 
This can be seen by expanding (2.14), 

5F - - dF OF 
°=^-« + 9G' G + ^' r * + "" <2 ' 15) 



with all the derivatives evaluated at = 0o, G — 0, and r# = 0. There will be 

SF 



a solution for as long as the operator 4? in (2.15) is an isomorphism between 



★ In the appendix we sketch a finite dimensional version of the theorem, which illustrates the 
relevant points. 
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two Banach spaces H\ and H2, and the two functions ^ and belong to the 
space H2. The choice of particular function spaces depends on the system under 
consideration. However, roughly speaking, we can see that this will be true in 
general given that the flat space solution 0o is a minimum of the energy functional 
(2.12) , which is equivalent to 

8F 

— -5<j)>0. (2.16) 

d( P (0o,O,O) 

Hence ^ has no zero modes and is invertable. In the next section we indicate how 
to choose appropriate function spaces for global monopoles. 

The OV equation implied that (pf + p) oc at a horizon. Above, we found 
that this same condition was needed to integrate out the metric coefficients A 
and B from the action. This allowed us to use the existence of non-gravitating 
solutions to imply via the implicit function theorem the existence of gravitating 
solutions and solutions with horizons. If we take a theory, such as Q-balls, in 
which, as we will see below, A and B cannot be eliminated from the action, then 
to use the implicit function theorem, one would have to compute the variation 
including all the dependent functions, <p, A and B. But knowledge of the flat space 
solutions gives us no information analogous to (2.16) about variations in the A or 
B directions, so the argument can't proceed. 



3. Global Monopoles 

In this section we demonstrate the use of the implicit function theorem and 
selection of appropriate function spaces for global monopoles. The matter field 
theory for the basic global monopole is given by an 5*0(3) invariant Lagrangian 
for a triplet of scalar fields a , 

C = ±V> a V M a - ±A (0 ffl a - v 2 f , (3.1) 

where V p is the covariant derivative operator. The scalar field configuration for 
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the monopole has the spherically symmetric form 



<f = V (p(r)f a . (3.2) 

For solutions without horizons <f>(r) interpolates between at the origin and 1 
at infinity. Evaluated on such field configurations (with the covariant derivative 
operator appropriate for the spherically symmetric metric (2.1)) the lagrangian has 
the form L m = j^K + U, where the kinetic and potential terms are given by 



K = U=^- + V(0 2 - I) 2 , (3.3) 



Here = dip /dr. The equations of motion for the metric coefficients are 

m'(r) = Anr 2 (\K + U), = WirGrK. (3.4) 

A \AB) 

The flat space global monopole solution has the following asymptotic behavior 

r ar, r — > 0; 

Mr) ~ r ! (3.5) 

where a and b are constants (the slope a at the origin must be determined numer- 
ically). From (3.5) and (3.3), one can see that the energy density for the global 
monopole falls off only as 1/r 2 , so that the total energy of a global monopole 
diverges, 

lim m(r) = 4Trv 2 r. (3.6) 

r— >oo 

Hence the spacetime of a global monopole is not asymptotic to flat spacetime, but 
rather to flat spacetime minus a missing solid angle [16], 

lim \ = 1 - 8nGv 2 . (3.7) 

In order to avoid a horizon at large radius (which is not of the sort we are interested 
in), we will keep 8itGv 2 < 1. 
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The quantity ^ in (2.15) for the global monopole is given by 

It 1 * = ~Tr ( r2 Tr S *) + < 2 + ^ W " 2 J> «* (3 ' 8) 

Here we have rescale lengths by a factor V At> 2 . The variations ^ and evaluated 
on the background solution can be seen to have the forms 

f r, r — > 0; ^ (const, r — ► 0; 

(3.9) 



<9G t ^, r -> oo, 9r ft \ ^, r -> oo. 

If we take the variation 50 to have the assymptotic behavior 

( const, r — >■ 0; 
50- <^ , (3.10) 
If, r^oo, 

(with the standard I? norm in three dimensions), then we can accomodate the 
variations induced by (3.9). This can be seen by examining the asymptotic behavior 
of |? in (3.8). We then have to show that the operator L = j+ is an isomorphism 
between these spaces. Since the operator is elliptic, this will be the case if neither 
it nor its adjoint have zero modes. Suppose that L has a zero mode, then we can 
write 

00 



Integration by parts yields 



= -r 2 f±f 
dr 



d „ x0 5 2 U 



oo 






f drr 2 | 







The boundary term vanishes for functions / having the behavior (3.10). Equation 
(3.12) then leads to a contradiction if 

w - -° (3 - 13) 

d( P (00,0,0) 

holds everywhere. We have checked numerically that (3.13) is satisfied for the flat 
space monopole. Therefore the operator L, which is self-adjoint has no zero-modes. 
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4. Examples 



Three examples of field configurations which allow horizons inside are SkyrmionsJ 
gauge monopoles, and global monopoles. These three examples span a range of 
types: gauge monopoles have both a long range magnetic field and topological 
winding, global monopoles have only the topological constraint, and the Skyrmion 
field winds but is not topological. These all have L m (<f)) of the form (2.11), and 
so satisfy the condition (p + Pf)\ rH = at a horizon. The implicit function the- 
orem argument shows that solutions with hair exist for G and rjj in some range 
about zero, but gives no information about how large this range is. One can de- 
duce more information about the range from arguments based on the traditional 
positive 'no-hair' integrals, which we do below in Section 5. 

Field configurations which cannot support horizons include Q-Balls [1] and 
boson stars [2]. Q-Balls are star type configurations that exist without gravity [1], 
but, as we will see, fail to satisy the condition (p + Pf)\ rH = at a horizon. The 
simplest Q-balls occur in the theory of a single complex scalar field [1] . The Q-ball 
field has the form = f(r)e~ lujt where f(r) vanishes at infinity. The lagrangian 
evaluated on such configurations is 

L Q = ^(f f ) 2 + l(m 2 -^f) + U(n (4.1) 

where the mass-term in the potential has been separated out. The frequency uj 
must satisfy oo 2 > m 2 for stability. From the definition of the stress tensor we then 
have 

2 SL m 2 5L m 1 l 2 1 2 1 2 (a o\ 

We see that to satisfy {p + Pf)\ rH — 0, / must vanish at a horizon*. But this 
means that the field is in its vacuum both at the horizon and at infinity, which is 
not a Q-Ball type solution. 



* We assume that the volume element \[AB is well behaved at a horizon, which implies that 
B ~ r — th near the horizon. 
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Boson stars (see [12] for a review) are localised scalar field configurations which 
exist only with gravity. The matter lagrangian again has the form (4.2) (with differ- 
ent potential terms and with to 2 > m 2 ). Hence Boson stars satisfy (p + Pt)\ rH = 
only for /(r#) = 0, implying again that the field be in its vacuum at the horizon, 
as well as at infinity. 

A third example which probably does not allow hair is the Abelian-Higgs model 
[17]. If the scalar field has the form f(r) and the gauge field is given by A t (r), 
then the matter lagrangian is 

Lah = ^(/') 2 " ^K) 2 - ^<W/ 2 + \{f - v 2 ) 2 (4.3) 

This again is not of the form (2.11), and satisfying (p + p?)\ r = requires that 
A\f 2 = at r = r#. While this in itself is not enough to rule out solutions, it 
clearly makes it "harder" to satisfy the equations given this additional condition 
on the fields. Indeed, the 'no-hair' integrals discussed in section 5 further imply 
that if At(rjj) = 0, then the fields are in their vacuum states everywhere outside 
the horizon. Adler and Pearson [17] explicitly analyzed the Einstein equation for 
this system further, and have shown that this is indeed the case. 

Finally, it is interesting to think about the case of a Coulombic electric field 
due to a point charge. This is outside the framework of the present discussion, 
because the non-gravitating configurations are singular, At = q/r. However, the 
Reissner-Nordstrom charged black holes are solutions with nonzero, nonvacuum, 
regular matter fields outside the horizon^. In this case, it is easy to check that the 
E&M Lagrangian reduces to 

Lem = -^{A' t ) 2 (4.4) 

which has the form (2.11) and that, in fact, the combination Pf + p vanishes 
everywhere. 

f Visser [18] has independently studied the condition (p +Pf)\ rH = in the context of various 
recent black hole solutions in field theories, such as dilatons and axions, coupled to gravity. 
He has also looked at the thermodynamics of such solutions. 
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In looking at these various examples, one notices that different kinds of mass 
terms play quite different roles. A "true" mass, or any potential U which is inde- 
pendent of the metric, makes no contribution to the sum pf + p, as in Inflation. A 
dynamical mass which comes from the coupling to the time component of a gauge 
potential, contributes a term to p? + p oc ^f 2 A^, which tends to rule out hair. A 
dynamical mass which comes from coupling to the spatial components of a gauge 
field contributes zero, and contributes a winding term oc ^ to p^ — p?, which is 
important in the OV equation (2.4). 



5. 'No-Hair' Integrals 



It is interesting to see how the black hole solutions discussed above avoid being 
ruled out by standard 'no-hair' arguments. In the case of extended gauge monopole 
solutions, this was discussed in ref. [7]. We will see that Skyrmions and global 
monopoles escape in basically the same way. Necessary conditions for the existence 
of black hole solutions in a given field theory can be derived by constructing energy 
integrals from the equations of motion (see e.g. [17,19]). If the action in the region 
outside the horizon is given by 



oo 

- J drJ(r), 



(5.1) 



an extremum occurs when 



(5.2) 



d SJ _5J 

dr 5(f)' 50' 

with the boundary conditions SJ/S(p' = at r = r# and the fields going to their 
vacuum values at infinity. Therefore 



dr 



5J_ 

5<j)' 



+ (d>- 



{<P - (Poo) 



SJ_ 







(5.3) 



Consider the case at hand (2.12), where 5* = — E and J is the positive definite 
integrand. Since we are assuming that regular solutions exist when G = rjj = 0, 
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the above is true with = and e y = 1 in J. Since typically the gradient term 
in the integrand is of the form C 2 (0)(0 / ) 2 , this requires that as r ranges from zero 
to oo, there are positive and negative contributions to the potential (the second) 
term in the integrand. Now, if the lower limit is taken to be rjj, there is still a 
possibility for positive and negative contributions to sum to zero above, if rjj is 
small enough. This point was discussed in [7] in reference to gauge monopoles, 
noting that the fields had to be Reissner- Nordstrom outside the horizon if rjj were 
sufficiently large. For Skyrmions, the structure of the no-hair integrals depends 
on what the response is of the Skyrmion field to gravity. But assuming that the 
effect of gravity is to further concentrate the energy density, again there will be 
a critical value of rjj, such that if the horizon is larger, the field must be in its 
vacuum outside the horizon. On the other hand, global monopoles have no such 
restriction on the value of rjj. 

Acknowledgements: We would like to thank Karen Uhlenbeck for helpful and in- 
formative discussions and the Aspen Center for Physics for its hospitality during 
part of this work. 

APPENDIX A 

Here we recall the arguement for the implicit function theorem for a system of 
N equations in N unknowns, and the limit as N becomes a continuous variable. 
Let g be the independent variable, and ir ,i — 1, ...,N be N dependent variables. 
(These are numbers, not functions.) We seek solutions ir = <j>i(g) to the system 

F j (7r,g)=0,j = l,...,N , (Al) 

given that <pi is a solution when g — 0, F((f>i , 0) = 0. Let ir — <f)i = 5ir and 
denote the matrix of first derivatives with respect to the independent variables by 
Oji = —-frf, evaluated at 4>io,g = 0. Then Taylor expanding the equation F — 0, 
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to linear order one needs to solve 



f)F- 

Jt 5n = -^--9 (A.2) 

There is a solution 5tc for any "source" on the right hand side of (A.2) if the matrix 
Oij has no zero eigenvectors, i.e., 

OijvV + , for all v i (A.3) 



For an implicit functional theorem, we would like the limit where the discrete 
index i becomes a continuous variable x, with Fj — > F(x), n — > (f>(x). Let {Pi(x)} 
be a set of basis functions, and let 4>(x) = SjAjPj(x) and 5(/>(x) = T<i5AiPi(x). 
Then in this limit, 

E -f h - 1 d * TM sm = E, S M - (A ' 4) 

Hence for a solution one needs that this last quantity, evaluated at the known 
solution, has no zero modes. In the main part of the paper, this condition was 
met since the second variation of the energy functional was nonzero, at the non- 
gravitating solutions. 
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